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It has been claimed that an impurity hopping between two sites in a metal can lead to non- 
Fermi liquid behavior if the bias between the two sites is tuned to zero. Here it is shown that 
several extra tunnelling processes are allowed; those involving bare impurity tunnelling as well as 
two-electron-assisted tunnelling play important roles. These terms generally drive the system away 
from the intermediate coupling non-Fermi-liquid fixed point, which has the character of the two 
channel Kondo problem. In order to find the non-Fermi liquid behavior a combination of these 
terms must be small over a range of energy scales. The prospects for obtaining these conditions 
are studied by analyzing both the weak coupling and the intermediate coupling behavior, together 
with the connecting crossover region. Unfortunately, the conclusions are rather pessimistic. The 
phase diagram and the important crossover are analyzed via bosonization and refermionization and 
some of the results compared to those obtained from conformal field theory. Finally, an additional 
non-Fermi liquid fixed point is found, a multicritical point with two relevant directions even in the 
unbiased case. 



I. INTRODUCTION 

A two-channel Kondo impurity in a metal is, per- 
haps, the most promising impurity system to exhibit low 
temperature non-Fermi-liquid behavior experimentally. 
This hope rests on the relatively few degrees of freedom 
involved (a local spin doublet coupled to two degenerate 
"flavors" of spin-^ fermions) with the symmetry group 
under which they transform being not very complicated 
[SU(2) fi avor x SU(2) spin x U(l) charge \. As a result one 
might hope that it would not be difficult to find a system 
whose low temperature behavior would be described by 
the two-channel Kondo system. 

Yet, for more than a decade since it was first intro- 
duced by Nozieres and Blandiro, no experimental re- 
alization of this model has been conclusively demon- 
strated. The difficulty lies in the fact that the non- 
Fermi-liquid fixed point is unstable to various symmetry- 
breaking processes which turn out to be present in real 
experimental situations. For example, Nozieres and 
Blandintj pointed out that anisotropy between the two 
flavor channels, caused by lattice effects, would destroy 
the non-Fermi-liquid ground state. Since this anisotropy 
could not be made to vanish for the known cases, the 
search for non-Fermi-liquid behavior with, conventional 
spin-Kondo systems was suspended. CoxB pointed out, 
however, that under certain symmetry conditions, local 
quadrupolar degrees of freedom could result in two chan- 
nel Kondo-likc coupling; unfortunately, dilute impurity 
systems of this type have provenJiard to make. 

Later, Vladar and ZawadowskH suggested that a non- 
magnetic impurity tunnelling between two sites in a 
metal could be modeled as a two-channel Kondo sys- 
tem in which the roles of the channels and the spins in 
the original formulation are interchanged. In this sys- 
tem the spin of the electron plays the role of the "flavor 
channels" , so that the anisotropy between "channels" is 



no longer an issue since in zero external magnetic field 
the spin-up and the spin-down electrons are degener- 
ate. This led Vladar and ZawadowskiQ to predict non- 
Fermipiliquid behavior in such a system. Recently Ralph 
et alau have interpreted low temperature tunnelling data 
in very small metallic contacts in terms of two-channel 
Kondo-like physics. Their measurements are claimed to 
be consistent with certain exact results obtained by Af- 
fleck and LudwigdjJput at this point, the interpretation 
is still controversial!! . 

Unfortunately, the mapping of the two-site impurity 
to the two-channel Kondo (2CK) system is far from 
exact, and, even with no anisotropy between "chan- 
nels" (the spin-up and spin-down electrons) there are 
other processes present in the tunnelling impurity sys- 
tem which are relevant and hence generically destroy the 
non-Fermi-liquid ground state. These processes, which 
cannot be neglected have not been treated adequately in 
the literatureotm. 

In this paper we carefully consider the mapping be- 
tween the tunnelling impurity system and the two- 
channel Kondo (2CK) problem. We first analyze the 
behavior when the impurity tunnelling starts to become 
important as the temperature is lowered; this is analo- 
gous to the weak coupling regime of the Kondo problem 
and is needed in order to understand which tunnelling 
processes will dominate at low temperatures. We then 
analyze the intermediate coupling behavior to investigate 
whether there are parameter regimes in which the system 
will be governed by the 2CK fixed point over a reasonable 
range of temperatures - as would be needed to observe 
the non-Fermi liquid behavior experimentally. 

Unfortunately, even in the optimal case of an impurity 
tunnelling between two identical sites so that the sys- 
tem has an extra Zi symmetry, we find that the 2CK 
fixed point is only accessible if two tunnelling processes, 
which are very different physically, nearly exactly cancel. 
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Generically, the system will exhibit Fermi liquid behavior 
at low temperatures. 

We analyze the behavior near the 2CK point in detail, 
focusing on the connection between the physical opera- 
tors and those that appear as "natural" operators in the 
2CK language. It is found that, on the critical manifold 
which can be obtained in the symmetric impurity prob- 
lem by adjusting one parameter, there are four leading 
irrelevant directions in contrast to the behavior for the 
pure 2CK problem analyzed by Sengupta and Georgestj. 
Various symmetry breaking terms are also studied and 
our results recover those derived by Affleck et atliB with 
conformal field theory techniques. 

A somewhat surprising feature emerges in the fuller 
phase diagram of the symmetric impurity model: a sec- 
ond fixed point which exhibits non- Fermi liquid behavior, 
albeit one with two relevant directions in the Zi symmet- 
ric case. 



A. Outline 

In the remainder of the Introduction we motivate the 
form of the Hamiltonian in which we focus and inter- 
pret the various terms that should appear. In Section II 
we derive an effective Hamiltonian that we will study and 
analyze its symmetries, while in Section III the weak cou- 
pling analysis is outlined. In Section IV the behavior near 
the intermediate coupling 2CK fixed point and, for com- 
pleteness, the various symmetry breaking operators near 
the 2CK fixed point are studied. In Section V we discuss 
the existence of an extra novel fixed point at intermedi- 
ate coupling. In Section VI, we discuss the accessibility 
of the 2CK fixed point and draw our conclusions. Finally, 
in the Appendices, the details of the weak-coupling (Ap- 
pendix A) and the intermediate-coupling analysis (Ap- 
pendix B) are presented; the comparison of our results 
to those obtained by conformal field theory is made in 
Appendix C. 



B. Physical Picture 

The system we wish to describe is an impurity or heavy 
particle which can hop back and forth between two sites 
coupled to a bath of electrons. The two sites may or 
may not be equivalent but we will primarily focus on the 
symmetric (equivalent) case. The asymmetric case can 
readily be treated in a similar manner. The effects of 
the interaction of this impurity with the electrons can 
be manifested in a number of ways. First, the electrons 
will tend to screen the charge of the impurity. Thus the 
impurity will hop between the two sites carrying with it 
tightly bound electrons which can move fast enough to 
adjust to the position of the impurity; it is convenient 
to consider these to simply be part of the "dressed" im- 
purity particle. However, in addition, as the impurity 



moves it may also redistribute the low energy electronic 
excitations near the Fermi surface. Since we are inter- 
ested in the low energy physics, these processes must be 
treated directly. For simplicity we consider only s-wave 
scattering off the impurity. 

If an impurity of charge Ze hops between two well sepa- 
rated sites, then the Friedel sum rule relates the (s-wave) 
scattering phase shift off a static impurity, 6, to the elec- 
tronic charge that will be moved to screen the impurity 
as it moves adiabatically from one site to the other, via 
Z = 2(5/7r, with the factor of two due to the two spin 
species. Conversely, if the two sites are close together, 
one can still usefully speak of an effective charge Q (per 
spin) which plays an analogous role to Q — S/ir in the 
well separated case, but is no longer simply related ei- 
ther to scattering phase shifts or to the impurity charge. 
It will instead turn out to be exactly the "orthogonality 
catastrophe" exponent that determines the system size 
dependence of the overlap between the electronic ground 
states with the impurity at the two sites.E3 For simplicity, 
we will focus on Q in the range < Q < 1, correspond- 
ing to repulsive interactions. (As shown in reference |l4| , 
other ranges of the effective charge can be reduced to this 
case via a set of more complicated combined impurity- 
electron processes related to those we consider here.) 

The important processes, in addition to hopping of the 
(dressed) impurity by itself, will be those in which one 
or two low energy electrons move in the opposite direc- 
tion to that the impurity hops. These processes can, for 
Q > 1/2, reduce the effective charge that must relax to 
the new impurity position, thereby decreasing the orthog- 
onality between the pre- and post-hop configuration; this 
results in a larger amplitude for the combined impurity- 
electron process at low temperatures, relative to the sim- 
ple impurity hop process. [Note that in reference a 
sign error in the definition of Q resulted in the incorrect 
interpretation being given to these processes; this error 
does not affect the conclusions, just the interpretation of 
the combined hopping processes] . 

In order to proceed with the analysis of the low tem- 
perature behavior of interest it would appear to be im- 
portant to assess the relative magnitudes of the ampli- 
tudes, at some high energy scale, of the processes dis- 
cussed above. However, it has been claimedEZl that this 
is essentially impossible for strong electKon-impurity in- 
teractions. Indeed Kagan and Prokof 'evt3 have claimed 
that a sensible Hamiltonian cannot be written in terms 
of a simple two level system since the high energy elec- 
tronic degrees of freedom cannot be properly taken into 
account. Although there are indeed real difficulties here, 
it should nevertheless be possible to introduce effective 
amplitudes at some intermediate energy scale and then 
analyze the behavior of the system in the phase space 
of these effective parameters. In general, unless there 
are specific reasons to prevent it, one would expect that 
most combinations of parameters could, in principle, be 
realized. We thus approach the problem via this route 
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and start with an effective Hamiltonian at an intermedi- 
ate energy scaie, at, say, some fraction of the conduction 
bandwidth. As we shall see, the extra hopping processes 
will in any case be generated at low energies. 



II. EFFECTIVE HAMILTONIAN 

The important electronic degrees of freedom at low 
energies are those that interact with the impurity in one 
of its two positions "1" and "2" . These are just the s- 
wave conduction electrons around the positions of the 
two sites. Thus at each energy, there will be two im- 
portant electronic degrees of freedom per spin. However, 
unless the two sites are very far apart (in which case the 
impurity tunnelling rates will be negligible and hence not 
of interest), the two sets of s-wave electrons will not be 
orthogonal; this will play an important role in what fol- 
lows. If we label the two sets of s-wave electrons by 
their energy, e, measured from the Fermi surface, then 
for each e there is an (essentially) unique pair of linear 
combinations of the two s-wave states, that are orthonor- 
mal and transform into each other under interchange of 
the two sites. We label the annihilation operators of this 
orthonormal pair c\ e and C2 £ with anticommutation rela- 
tions 



|c ie ,cj e ,j = 2nS (e - e') S tj 



(1) 



with the "1" and "2" denoting the sites near which the 
wavefunction is larger. The impurity interacts with, sim- 
ply, the operators 



de 

Cl,2 = / 7^Ci j2e , 
Z7T 



(2) 



although in each position the impurity will couple to 
both c\ and c 2 due to the non-orthogonality of the two 
electronic s-wave states. With time reversal invariance, 
which we assume henceforth, the most general interaction 
with the impurity becomes (ignoring spin for now) 



Note the appearance of an effective electronic hopping 
term V2, caused by the scattering by the impurity; this 
will vanish if the sites are far apart, but in general will 
be comparable to the other terms. 

The first term in Eq(^) which is the average of the 
interaction over the two impurity positions, merely pro- 
duces a constant phase shift for both "1" and "2" elec- 
trons at the Fermi level and, combined with other oper- 
ators, gives rise only to irrelevant terms. Therefore we 
will ignore it at this point although in Section IV an ir- 
relevant operator it gives rise to will play a role in our 
discussion of the intermediate coupling behavior. 

With the effective hopping charge Q in the range 
[0,1], there-. are three hopping processes that must be 
considered!! 4 ]: 



Hhop 



dXdi 



An 



A i ( t t 

f ~T v c it C2 t + c n c n 

-A 2 c{ T C2 T ci i C 2 X 



(5) 



representing hopping of the (dressed) impurity, jointly 
with, respectively, 0, 1 and 2 electrons moving the oppo- 
site way. Although we might start at a high energy scale 
with negligible Ai and A2, these will be generated under 
renormalization and hence must be included. 

In order to analyze the renormalization group flows, 
it is convenient to approximate the conduction band s- 
wave electrons c\ t and C2 e by a linear dispersion with a 
cutoff, at short times, r c , roughly the inverse bandwidth. 
Then the interactions with the impurity will essentially 
be replaced by corresponding phase shifts, specifically V3 
replaced by an effective phase shift that we denote wQo; 
Qo will have the interpretation of an effective "charge" . 
We then have, after inserting powers of t c to make the 
couplings dimensionless and factors of tt for convenience, 



H — Ho + Hint + Hhop 



with 



Ci^Clo-e + 4ere C 2<xe 



(6) 
(7) 



U = d\d x 



+d\d 2 



(oS 
Vi (c 2 c 2 



c\c 2 



c\ci 



+ v 2 


[c\c 2 


+ c\c\ 


+ v 3 


( c i c i 


- c|c 2 


+ v 2 




+ c\c 2 


+ v 3 


(4c 2 


- c\ci 



(3) 



where are the annihilation operator of the impurity 
at the sites. Using the obvious identity d\d\ + d\d 2 = 1, 
Eq(|J) can be written as 



U = Vi (c\ci + 4c 2 ) + V 2 (c\c 2 + 4ci 
+ V3 (d\di - d\d 2 ^j {c\ci - 4c 2 



(4) 



and 



nQ [d\d 2 - d\dij (4 ( t C 2<t - 4<r C lc 

a 



(8) 



Hhop — d\d\ 



A , Ai ^ f 



2lTT c 2 

a 

+ A 2 2tTT c c\ 1 C 2 ^c\ |C2| 



h.c. 



(9) 



We have rescaled the electronic term V 2 to a coefficient 
y which will play the role of a "fugacity" for electronic 
hops. 
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A. Symmetries 

It is important at this stage to examine the symme- 
tries of Eq(||). In addition to time reversal, conservation 
of the electrons (c — + e 4 *c) , conservation of the impurity 
(d— > e^d), and SU{2) spin symmetry, the only other 
symmetry is interchange of the two sites and the corre- 
sponding electronic states (1 <-*■ 2). Note, however, that 
if the only hopping term had been Ai, and if y vanished, 
there would be an extra artificial symmetry d\ — > e l ^di, 
ci — > e 4 ^ci, di — > d2, C2 — > C2 corresponding to conser- 
vation of iVi = djdi + n ci and similarly N2 separately, 
where n Cl is the number of the "one" electrons which, in 
the absence of the channel mixing term y, are indepen- 
dent of the "two" electrons. As shown in reference 
even if the electronic states had been optimally chosen so 
that there was no mixing of "one" and "two" electrons 
at the Fermi energy, the energy dependence of scatter- 
ing off the impurity would generate extra mixing terms 
in Hint, that cannot simply be expressed in terms of c\ 
and C2- These would break the artificial symmetry and 
under renormalization generate a y (c\c2 + h.c^j mixing 

term even in the absence of impurity motion. Thus it 
is best to include y from the beginning. (The neglected 
energy dependent scattering terms will then not play an 
important role). In order to understand the difficulties 
of reaching the 2CK-like fixed point, this step is crucial. 

The artificial symmetry in the absence of the chan- 
nel mixing and Ao, A2 terms, corresponds to a con- 
served pseudo-spin N± — N2 which is the sum of the 
"z-components" of the impurity pseudo-spin d\d2 — d\d\ 



and an electronic pseudo-spin n c 



This pseudo-spin 



can play the role of spin for the two-channel Kondo effect 
and, indeed, under renormalization the system will fl 
to this intermediate coupling 2CK fixed point if Q > 
Unfortunately, there is no natural small parameter which, 
keeps the pseudo-spin symmetry breaking terms small.L£l 



B. Bosonization 

In order to carry out the renormalization group analy- 
sis for small bare impurity hopping rates, it is useful, as 
is standard, to bosonize the electronic degrees of free- 
dom, treating the electronic states c\,it as those of a 
one-dimensional system with two sets of right moving 
electrons with "wavevectors" Vp{k — kp) oc e. It is sim- 
plest to set the Fermi velocity, Vp = 1, and treat e like a 
wavevector index, defining 



1 



ere 



so that 



ct (x)c- (x)= 1 d ^ Ax) 



(10) 



(11) 



with j = 1,2 a =t,l and <&j a being the correspond- 
ing bosonic degrees of freedom.-where we have followed 



Emery and Kivelson's notationBS Only . 



(x = 0) 



couples to the impurity. Note that in the standard ex- 
pression Eq(pd]) the left hand side is normal ordered and 
therefore the (infinite) uniform charge density does not 
appear. Also corrections that vanish as r c — > are ne- 
glected; we will be careful to include the effects of extra 
terms when they play an important role. 

Since we will later need to be careful to have the proper 
anticommutation relations, we must insert extra factors 
of the form exp (iirN^), with = J dx** (x) M/ M (x), 
into some of the bosonized expressions to ensure anti- 
commutations of the different Fermi fields. These will 
not play a role as long as no spin- flip processes occur, 
and, for the time being we ignore them; the needed mod- 
ifications are spelled out in Appendix B. 

It is useful to define even and odd components of the 
Bose fields 
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the Hamiltonian then becomes 



V2 



dx 



I COS $ c 



(12) 



(13) 



A 

2lTT c 

A 



Ai 

47T7V 



E 



+ exp 



[iV2$ c 



h.c 



+ (cr+ exp [j-v/2 ($oT + $ol] 



where in all the coupling terms the bosonic fields are 
evaluated at x = and we use the impurity pseudo-spin 
operators 



and 



<j z = d\d2 — d\d\ 



<J+ = d\d\ , er_ = d\d2- 



(14) 



(15) 



The conduction electron part of the Hamiltonian can 
be written in terms of boson creation and annihilation 
operators <j)^ (e) with canonical commutation relation 



(6),0t ( e ')j =27^(6-0 



(16) 



via 



*« (x) 



and 



C°° df 



(17) 



(18) 



4 



which involves positive energy parts only. Here r c 1 is the 
energy cutoff and fi represents the various Bose fields, i.e. 
for Eq|[!), fi = (e t, e 4, o T, o J.). We see that Eq(g|) 
does not include any terms with or $ e |. Thus, up 
to operators that are irrelevant for weak coupling, the 
impurity is decoupled from the even boson fields. ■— ■ 

It is convenient, following Emery and KivelsorJIS, to 
decompose the odd field $ Q which couples to the impu- 
rity, into a spin and charge part, by 



dA (I 



V2 



($ c + a$ s ) 



(19) 



with a = ± for spin f, J., respectively. The second 
term in Eq(|l3|) becomes simply Qo<J z ^ (0). This term, 
which represents the difference in phase shifts for electron 
scattering off the two positions of the impurity, can be 
shifted away by a unitary transformation which changes 
the naive weak coupling scaling of the, bopping terms; 
this is the conventional approach usedBE3o'E£l to derive 
the weak coupling flows discussed in the next section. 

Although the even fields will not play much role for 
the time being, for later purposes we also introduce even 
fields $ ec , $ es by 



cr$ P 



(20) 



Note that any sum or difference of any two of the &j a , 
i.e. those that appear from operators bilinear in electron 
operators, can be written as a sum or difference of two of 
the fields <& s , $ c , $ e s, ®ec with coefficients roi unity; this 
enables the method of Emery and Kivelsonll3 to work. 



III. WEAK HOPPING ANALYSIS 

In order to connect the various amplitudes at the rel- 
atively high energy scale of the effective Hamiltonian 
(Eq(|l3|)) to their renormalized values at low energies we 
must analyze the weak coupling renormalization group 
(RG) flow equations for the amplitudes in Eq([l3|). The 
magnitudes of the various terms at the crossover scale 
to intermediate coupling will determine which regions of 
the initial parameter space can flow near to the 2CK fixed 
point. 

Following the procedure in reference |l4| we transform 
7i to UHW using the unitary operator 



U = cxp [-ia z Q <S> c 



(21) 



Subsequently we follow the RG approach described there 
and obtain the following flow equations for the various 
amplitudes, where for later convenience we introduce 



Q = ^ - Qo, 



(22) 



d} v ,2q- 2q 2 ) A + yAi + O (A J ) (23) 

^=Q-2g-2 9 2 )A 2 + y A 1 + 0(A^ 
dAi fl 



= - - 2<f Ax + 2y (A + A 2 ) + O (A 3 ) 



dl 



Ai (A + A 2 ) 



The important cross-terms in the first three equations 
in Eq(p3|) that are proportional to the electronic mixing 
term y, have a simple physical interpretation: they rep- 
resent the effects of an impurity and an electronic hop 
both occuring within a short time interval so that, at 
lower energies, this appears as simply the corresponding 
combined process. Note that we have not included O (y 2 ) 
terms in the above equations; the definition of these will 
depend on the RG procedure, and they will not qualita- 
tively change the behavior. Thus, in the spirit of focusing 
on the important processes and terms, we ignore therrO. 

Noting that q and y are constant to order O (A 2 ), to 
analyze the flow for weak hopping, we can safely set them 
to their initial values, go and Vo- Then the first three 
equations can be diagonalized exactly; the details are 
discussed in Appendix A. The RG eigenvalues for the 
hopping terms, about the zero hopping fixed line are 



(24) 



A °=2 



which lies in the range (—1/2, 1/2): 



We now note that for < Qo < 1, corresponding to 
1 90 1 < 5' a t l eas t two eigenvalues are positive so that 
impurity hopping is always relevant (leading to the con- 
clusion of absence of impurity localizations) ; likewise for 
other ranges of Qo, there will always be at least two rele- 
vant hopping processes if there is only s-wave scattering 
off the impurity. 

The Kondo temperature, T^-, is the energy scale at 
which the first of the impurity hopping processes becomes 
of order unity- i.e. of order the renormalized bandwidth^ 

The system considered by Vladar and ZawadowskiH 
and Vladar et a£l, essentially amounts to neglecting A 2 
and y; for small \q\, which will turn out to be the most 
interesting case, this misses part of the physics. The 
reason is simply that they neglect one relevant operator 
A 2 which mixes with the other two to give the correct 
eigenvalues (Ecj24|). Furthermore, as will be seen later, 
non-vanishing values of A 2 and y are crucial to give the 
correct renormalization flows close to the intermediate 
coupling fixed point. 

If we only kept Q and Ai non-zero, their weak cou- 
pling flows would be (up to coefficients) like those for J z 
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and J± for the conventional Kondo problem. The Kondo 
scale is then simply 



Ai\- 

ir \w J 



(25) 



after reinserting factors of the bandwidth, W. (For the 

_ i 

special value Ai = 2Q, Tk ~ e like in the well known 
anti- ferromagnetic Heisenberg Kondo problem) . For this 
artificial case, at scales below Tk the novel two chan- 
nel Kondo physics will indeed appear, as we will discuss 
later. This results from the approach, at low energies, of 
the system to the intermediate coupling 2CK fixed point. 

Unfortunately, the breaking of the artificial pseudo- 
spin symmetry leads to the appearance of terms which 
generically drive the flow away from the 2CK fixed point. 
In order to analyze whether the system can get near to 
the 2CK fixed point — the prerequisite for observation of 
non-Fermi liquid behavior — , we must be able to identify 
the operators near the 2CK fixed point in terms of the 
original terms in the Hamiltonian; the magnitude of the 
operators, in particular the relevant ones, can then be de- 
termined, roughly, by "matching" the coefficients at the 
crossover scale, Tk, between the weak and intermediate 
coupling regimes. 

From Appendix A, we see that the crossover tempera- 
ture, Tk, will generally be a complicated function of the 
original parameters. Before examining the magnitude of 
the various important terms at Tk, we turn to the be- 
havior near the 2CK fixed point; this will tell us which 
terms need to be small at scales of order Tk for 2CK 
behavior to obtain. 



IV. INTERMEDIATE COUPLING ANALYSIS 
AND TWO CHANNEL KONDO FIXED POINT 

From the weak coupling flow equations in Eq(|23|) , it is 
apparent that, in the absence of electronic mixing (y=0) 
there is a special value of the effective impurity charge 
Q: for Q = \, corresponding to q = 0, Ao, Ai and A2 all 
scale in the same way with eigenvalue A = Ao = A± = i . 
By analogy to the Toulouse limit of the conventional one 
channel Kondo problem, it is thus natural to look for a 
solvable point that corresponds to y = q = 0, inspired 
by the observation that free Fermi fields have scaling di- 
mension of o and thus might be used to represent all 
of the hopping terms that appear for q = 0. This has 
been carried out recently by Emery and KivelsonEj who 
"refermionize" the bosonized operators that appear in 
the Hamiltonian enabling the computation of the scaling 
of the various important operators. 

By examining the weak coupling flows, it is apparent 
that special behavior might occur when Ao + A2 = 0. 
Defining 



we see that, at least to the order included in Eq(|2^), for 
y = and A + — 0, q flows to zero, y is not generated and 
one might hope that the flow would be towards the 2CK 
fixed point. Indeed, the intermediate coupling analysis 
shows that this can occur, even if 



A_ = A n - A 2 



(27) 



is non-zero so that the artificial pseudo-spin symmetry is 
broken. 

Physically the role of A + is very surprising. The pro- 
cesses represented by Ao and A2 are very different and 
the definitions which make them dimensionless are clearly 
cutoff dependent. Thus the special critical behavior must 
not in general occur exactly at A + = 0, since the loca- 
tion — but not the existence — of the critical manifold 
will be affected by irrelevant operators. In particular, 
from the weak coupling flows we can see that a non-zero 
q combined with A_ ^ will generate A + , thus a "bare" 
A + that is non-zero will be needed for the flow to go to 
the 2CK fixed point asymptotically. 

In order to understand the intermediate coupling be- 
havior and th£. special role of Q = \, following Emery 
and Kivelsor£3 and the analogous Toulouse limited of 
the conventional Kondo problem, we perform a unitary 
transformation with 



U = exp 



which transforms the Hamiltonian to Ti, = UTLU^ 

Ai A_ 
n = Ho + a x cos $ s + a y sin <I> C 



2nT r 
A, 



2jTT r 



(28) 



(29) 



, -o~ x cos $ c + 2y cos $ c cos 4? s qo~ z — — 

2lTT c OX 

ud$ ec w d$ ec 

--5 7T a x C0S ®c-Z— 

it ox 2ir ox 



where we have combined the terms 



2-KT c 



\<j+e " + er_e 



2itt c 



A_ 

27TT- 



A_ 



■a v sin $ c + x cos $ c , 

27TT C 



(30) 



A + = Ar 



A 2 , 



(26) 



and abbreviated <I> A1 (x — 0) simply by <i> M . We have 
also reintroduced some of the coupling terms to the 
even fields, in particular the marginal term u d g^ c that 
arises from the impurity-position independent part of 
the electron-impurity scattering {V\ in Eq(|])) and the 
term wa x cos $ c 9 q£ c , which arises from the combination 
of u and impurity hopping terms A + and is irrelevant for 
weak hopping, will play roles in our analysis. An addi- 
tional irrelevant term, a y sin $ s a J^ s couples the impurity 
to the electronic spin degrees of freedom, but does not 
feed back to the other operators and thus we ignore it 
for the time being. Its effect will be discussed further in 
Appendix C. 



G 



A. Symmetries 



4- Spin reversal 



Since we expect that the symmetries will play an im- 
portant role, we should examine what the original sym- 
metries correspond to in TL and ensure that there are 
no extra symmetries which might have arisen from the 
discarding of operators which were naively irrelevant, 
since, as shown in reference [l4|, such procedures, espe- 
cially when combined with "large" transformations, such 
as Eq(pq), can be dangerous. 



Flipping electron spins is simply <3? s 



and 



1. Gauge invariance and spin conservation 



5. Time reversa 



Since there are no spin-flip processes, separate gauge 
transformations can be made for each spin <frj a (x) — > 



^ja {%) + a , corresponding to $ e 



$ eff + V28 a . This 
does not play much role as the <I> ec -field enters only as a 
derivative ^J^P and the $ es -field decouples from the im- 
purity at the level at which we work (there is feedback, 
under renormalization from other operators involving $ e 
but these only modify pre-existing terms in TL). However, 
if the z-component of electron spin is not conserved, then 
only the symmetry $ ec — > 4> ec + 9 ec remains. 



2. Periodicity 

The definition of the Fermi fields, Eq(|l0[), implies that 
shifting any by 2ir should leave TL unchanged. De- 
pending on whether one or both spin components are 
shifted, this implies that (ignoring shifts in <E> ec ) 



(31) 



and 



$ s -> $ s + 2ir (32) 
are independent symmetries, as is the combination 

$ c -» $ c + 7T , $ s -> $ s + 7T (33) 



7T 



V 



'Vi 



with the necessity for the simultaneous transformation of 
o~ X y resulting from the unitary transformation, U, which 



involves $ r so that $ r 



7r introduces an extra 



exp (— ^cr 2 ) = — ia z factor into U yielding a x _ y — > —cr x ,y 
in TL. 



3. Interchange 

Interchanging sites one and two is equivalent to 

$ c -» -$ c , $ s -» (34) 



with $ e c,es unchanged. 



Time reversal transformations change ingoing to out- 
going waves, thereby yielding x — > —x, i — > —i, all 
^ (x) — ► — <I> M (— x) and a y — > — o^. Note that here we 
are not time reversing the spins. 



6. Artificial extra symmetries 



We now see that, indeed, TL in Eq(|29|) with all coeffi- 
cients non-zero, does not have any artificial symmetries. 
But as seen earlier, an artificial extra pseudo-spin sym- 
metry is possible: pseudo-spin conservation mod 2 corre- 
sponds to $ c — > <f> c + 7r, and more generally full pseudo- 
spin symmetry corresponds to independence of the "one" 
and "two" electrons, i.e. <& c — » $ c + 6* c with any 6* c . The 
terms w, A + , A_ and y all violate this, and it can read- 
ily be seen in the representation of TL of Eq(f29|) that y 
combined with Ai generates A+, and q combined with 
A_ generates A + , as expected. 

In the representation of Eq(^) we see that there is 
another possible artificial symmetry: If w, A+, y and q 
are all zero, then 



4\ 



7T - 3> 



(35) 



becomes a symmetry. This, as we shall see, restricts the 
system automatically to the stable critical manifold of 
the 2CK fixed point. But note that because of the uni- 
tary transformation of Eq(p8|) <& c — > 7r — $ c does not 
correspond to a realizable symmetry in terms of the orig- 
inal variables since it mixes hops involving the impurity 
alone and those involving the impurity together with two 
electrons. (Indeed many other irrelevant terms neglected 
in TL will also violate this artificial symmetry.) 
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B. Refermionization 

If A_|_, q, y and the other operators neglected in Ti all 
vanish, then the system will still exhibit the novel in- 
termediate coupling 2CK behavior. As we shall see, the 
only relevant operator near the 2CK fixed point, which 
is consistent with the true symmetries of the impurity 
hopping between two equivalent sites, is A + , thus only 
one combination of physical quantities needs to be ad- 
justed to obtain the 2CK behavior. Unfortunately, this 
is a combination which is not naturally small. 

The behavior near the 2CK fixed point,can most easily 
be found following Emery and Kivelsontij by "refermion- 
izing" the Bose fields $ c , $ s and $ ec that appear in H 
in Eq(pg|), noting that exp«$ M , (with <j>^ properly nor- 
malized) is like some pseudo-Fermi field. The details are 
discussed in Appendix B. 

Crudely, for /i = c, s, ec, es, each field e**f is replaced 
by a new Fermi field, W^, and cr„ by a local Fermi field 
d, with appropriate factors of e J7rAffl and e I7rJVd to give the 
correct anticommutation relations. The symmetries can 
most easily be seen, and the Hamiltonian simplified, by 
writing the new Fermi fields in terms of a set of Majorana 
(hermitian) fermions: 



d = 



1 



0) 



- (7 + ^) 
-7= {a» + i&n) ■ 



(36) 



Note that 



er 2 = 2i7<5. 



(37) 



The symmetry restrictions can now be examined in terms 
of these variables; the details are given in Appendix B. 
The periodicity of the Bose fields <i> c and <i> s simply im- 
plies that only terms with an even number of Fermi fields 
can appear in the Hamiltonian. Gauge invariance implies 
that can only appear as ^\ c ^ ec and the z-component 
of spin conservation that ^ es cannot appear in the ab- 
sence of magnetic fields in the x- or y-direction. Spin re- 
versal, because of the role of the ordering operators, takes 
^ s — > — and \I/es — * — ^Lj implying that a s and a es 
by themselves are excluded by spin reversal symmetry 
Interchange symmetry takes ^> c — ► — ^J, — ► — and 
d — > d) thereby requiring that a c and S must appear to- 
gether. Finally, time reversal takes x — > —x, i — > —i and 
$ — > — $, allowing only real coefficients of ^ operators, 
and hence forbidding terms like 17a. The Hamiltonian 
becomes 



TL = H.q 



(Ai 7 ft + A_ 7/ 3 C + A+Sa c ) (38) 



\/2TrT c 

- AiTy^5a c l3 s + 4irq~f5a c f3 c 

+ 2iua ec (3 ec + w\ / 2TrT c 5a c a e c{3ec 

with Ho the kinetic energy of the four (eight Majorana) 
new Fermi fields, ^ fl . With the full symmetries of the 



system, the other five fields (a s and the ec-, es-fields), 
cannot appear in the couplings to the impurity, except 
in relatively innocuous forms, involving the simple po- 
tential coupling to the average position of the impurity, 
ia ec (3 ec and combinations of this with other terms, as well 
as the irrelevant term 5a s a es (3 es , which will be discussed 
in Appendix C. Note that other potentially important op- 
erators, like 8(3 s a c j3 c and iS dot g^ are excluded by time 
reversal invariance. 

The original 2CK problem studied— iby Emery and 
Kivelsonlla and Sengupta and Georgest2l corresponds to 
A_ = A + = y = w = 0. In our case non-zero A_ can 
be important by observing that "half" of the impurity, 
7, couples to both {3 C and /3 S , thus it is convenient to 
rediagonalize and make linear combinations of these, pi 
and fix yielding, with 



A K = J&( + A* 



It = H.o 



(39) 



(40) 



+ -j== (AkiPi + A+Sa c ) 

+ AirqjSctcPx + 4,^yjSa c Pi 

+ 2iua ec p ec + V2iTT c w5a c a ec p e 



where y and q ar e linear combinations of the original y 
and q (see Eq(B7) in Appendix B for details). The above 
rediagonalization of /3 C and /3 S roughly corresponds to a 
rotation of "spin" axes in the conventional 2CK language. 

From the electronic kinetic energy Tio, the a's and (3's 
all scale, with time scale r, as r _ 2. If all the couplings 
are small, then the anti-commutation relations of 7 and 
S imply that they are dimensionless so that A + and Ak 
scale as r~2, while q and y are marginal as from the 
weak coupling analysis of Section II (Eqj23|)) and w is 
irrelevant. 



C. Two channel Kondo fixed point and flows 

When y = q = A + = w = 0, the Hamiltonian in 
Eq(JIo|) corresponds ±0 the Toulouse limit analyzed by 
Emery and Kivelsonta. As a free fermion system it can 
be analyzed straightforwardly. In this limit "half" of the 
impurity, 5, is uncoupled from the electrons and thus has 
no dynamics, while the other "half", 7, gets dynamics 
from coupling to the electrons with correlations at large 
imaginary times 



(T t7 (t) 7 (0)>~~. 

T 



(41) 



Together these yield the non-Fermi liquid 2CK behavior 



(T T a z (r) a z (0)) ~ -. 

T 



(42) 
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It is important to note that@ the solvable Hamiltonian 
is not generally at the 2CK fixed point. Indeed, the cor- 
relations are readily seen to exhibit crossover from weak 
coupling behavior, (<TzO~ z ) ~ const, to the non-Fermi liq- 
uid behavior of Eq(M2) for r > A k 2 t c . 

The 2CK fixed point, formally, corresponds to Ak — > 
oo . It is more convenient, however, to allow instead the 
normalization of 7 to change, corresponding to letting the 
coefficient of the J 7<9 r 7 dr in the Lagrangian vary. At 
the fixed point, this coefficient, say <? 7 , will be zero, while 
Ak becomes a constant; the correlations of 7 are then 
simply the inverse (in frequ ency space) of those of i.e. 
a pure power law Eq(BlC). To connect the two regimes 
together, one could choose to renormalize so that, for ex- 

A 2 

ample, - xg- + g 1 = 1 [a particularly convenient choice; 



dA 



+ 



dl 2 A ~ 



2yA 



wu 



K 



V2(l + u 2 )'- 



(43) 



(see Eq(B12)], by rescaling 7 under renormalization by 
a Ak -dependent amount. Details about this procedure 
are given in Appendix B. [Note, however, that the re- 
sulting fixed point Hamiltonian (with Ak = V47r) will 
not have the pseudo-spin SU (2) symmetry of the pure 
2CK problem. This is because the RG approach we have 
implemented here, including the unitary transformation 
of Eq(p8|), is inherently anisotropic in channels (equiva- 
lent to spin of conventional 2CK problem) . In Appendix 
C we will show that this anisotropy will not affect the 
results.] 

At the 2CK fixed point, the above scaling implies 
7 ~ t~5 i while S is still dimcnsionless so that the RG 
eigenvalues of the other operators can be read off im- 
mediately: + 1/2 for A_|_, the unique relevant operator 
consistent with the symmetries of the problem; -1/2 for 
q, y and w, which are the three leading irrelevant opera- 
tors discussed so far (the fourth will appear in Appendix 
C); and for u, which is marginal but redundant, in that 
it does not affect the impurity dynamics. 

The operator corresponding to q does not give rise to 
any terms which couple S linearly to the a's and /3's and 
hence will not generate A+. It is the leading irrelevant 
operator for .-the 2CK problem identified by Sengupta 
and GeorgesEH An artificial symmetry is responsible 
for its special role, indeed just the one discussed earlier: 
Eq(|35|), $ c — > 7r — $ c corresponds to the discrete sym- 
metry a c — > —a c and fix ~ > —fix which does not have a 
natural representation even in terms of $ c and $ s . But if 
present, this artificial symmetry excludes the generation 
of terms like A+, y and w, even in the presence of q. In 
fact, without these terms, the extra artificial symmetry 
is really an 0(2) symmetry in the {a c ,(3x} pair, consist- 
ing of a U(l) of rotations in the (a c , 0x) plane, combined 
(in a non-commutative way) with Zi , the usual site inter- 
change symmetry, which takes a c — > — a c and (3x — > Px- 
This is exactly analogous to the 0(2) symmetry present 
in the model treated, by Emery and Kivelsonlla and Sen- 
gupta and GeorgeaLO. 

In contrast, the irrelevant operators w and y, break 
the artificial symmetry and yield, at lowest o rder, the 
generation of the relevant operator A + (see Eq(B14)) 



consistent with expectations from weak coupling. This 
implies as stated earlier, that the critical point will not 
be exactly at A + =0. 

Of the three leading irrelevant operators, it should be 
noted that, although y has the same scaling dimension 
as q and w, it has a different role close to the 2CK fixed 
point. The reason is that it couples to the term ijPi, 
already present at the fixed point. But the Ak term at 
the fixed point suppresses fluctuations of /3j, causing the 
leading term in the correlations of /3j to vanish at the 
fixed point (with sub-dominant terms caused by g 1 ^ 0). 
As a result, unlike q and w which each yield. O (T\nT) 
contributions to the impurity specific heato — a key 
feature of a 2CK non-Fermi liquid — the singular part 
arising from y is only of O (T 3 In Tj , for temperatures 
T <C Tk . Thus only two of the above independent lead- 
ing irrelevant operators give leading singular specific heat 
corrections. In fact, there is a third one, involving only 
spin degrees of freedom, which is discussed in Appendix 
C. 

In Section VI, the conditions for accessibility of the 
2CK fixed point are analyzed. We turn here to further 
analysis of the behavior near the 2CK fixed point. 



D. Symmetry Breaking Operators 

Up to this point we have only dealt with an electron- 
impurity system which is invariant under Z2 (1 «-> 2 in- 
terchange) and spin SU (2) symmetry. However, in a re- 
alistic situation of an impurity in a metal, there will gen- 
erally be a non-zero, although possibly small asymmetry 
between the two sites which will break the Z2 symme- 
try. Furthermore, in the presence of a magnetic field, the 
equivalence between the two spin channels will be lost, 
leading to-a-situation similar to the anisotropic Kondo 
problemErcac 3 ]. As might be expected, in both cases, the 
symmetry breaking terms are relevant and in their pres- 
ence the system flows away from the 2CK fixed point. 
In this section, we will briefly comment on the effects of 
symmetry breaking terms close to the 2CK fixed point. 

It is clear from the discussion in the previous section 
that for an operator to be relevant close to the 2CK fixed 
point, it has to be of the form iSx where x is a Majorana 
fermion of scaling dimension 1/2. From the ten Majo- 
rana fermions (four pairs of a^, (3^ and 7, 6, all listed 
in Table Q) we can make nine such operators. Excluding 
iSa ec and i8(3 ec due to total electron number conserva- 
tion (which only allows a ec and (3 ec to appear together 
as a ec f3 ec ) we are left with seven possible terms. 

From the transformation properties under the discrete 
symmetries of the system, listed in Table Q, it can be seen 
that /3 C and (3 S have the same symmetries; indeed this is 
why the Ai and A_ terms in Eq(|38|) could be combined 







into the Ak term of Eq(ffO|). In the presence of small iS/3 c 
and iSP s terms, a small rotation of the (6, 7) pair as well 
as a small additional rotation of the (p c , (3 S ) pair can 
be performed to yield just a slightly modified iA^Pi 
term, and a single remaining relevant perturbation, the 
iSPx term. The extra operators (ijPx and id Pi) are thus 
"redundant" BJ Therefore at the 2CK fixed point, there 
are exactly six relevant operators, all with RG eigenvalue 
of 1/2, like A+. These, along with their symmetry prop- 
erties, are listed in Table ||. 

The first three operators in Tabic || do not break the 
spin SU(2) symmetry. Among these, the first, our fa- 
miliar iA + da c , is interchange and time-reversal symmet- 
ric, corresponding to the relevant part of the channel 
pseudo-spin operator S x . Correspondingly, the second, 
iSj, breaks interchange symmetry but is time reversal in- 
variant, corresponding to a S z operator. This will result 
from simple asymmetry between the impurity energies of 
the two sites, i.e. a o~ z term in the original Hamiltonian. 
Finally, i5(3 c (or, equivalently, i8/3 s ) breaks interchange 
and time reversal which makes it an imaginary operator, 
i.e. it is generated by a S y operator in the channel sector, 
corresponding to complex hopping matrix elements. 

The relevant spin SU(2) breaking operators are iSa s , 
iSa es and i8f3 es ; these correspond to joint electron- 
impurity hops accompanied by a spin flip or carrying 
electronic spin. They correspond to combinations of 

11C24,), o+ct.C2f, c+c||C2X and their her- 



C2T 



mit ian c onjugates and are discussed in Appendix B (see 
Eq flB16D and Eq (|B17|) ). The first corresponds to the 
"flavor" anisotropy .term in the conventional two chan- 



nel Kondo modetloa, and, being interchange and time 
reversal symmetric but odd under spin flip, is induced 
by a magnetic field in the z-direction. Interestingly, the 
remaining two relevant spin SU (2) breaking operators 
are odd under the Z2 interchange transformation. This 
means that in order to adjust them to zero in a non- 
zero external magnetic field, one would have to tunc also 
terms that break the interchange symmetry of the prob- 
lem, making any additional novel, finite magnetic field, 
non-Fermi liquid fixed points (analogous to that found in 
zero field for A + 3> Ai,A_), extremely hard to observe 
in a system that does not have interchange symmetry. 



V. ADDITIONAL INTERMEDIATE COUPLING 
FIXED POINT 

In the previous section we analyzed the behavior of the 
system close to the 2CK fixed point, that corresponds to 
the limit Ak ^> A+. There, we showed that 7, "half" 
of the impurity, acquired non-trivial dynamics, which es- 
sentially gave rise to the non-Fermi liquid behavior of the 
system. However, as is evident by examining Eq(|40|), it 
should be, in principle, possible to get the same type of 
non-Fermi behavior if the inequality above were reversed 
(A K < A+). 



Indeed, following the same arguments analyzed above, 
we see that when =q=y=w=u=0 and A + 7^ 
we have a critical point, which has an extra artificial U(l) 
symmetry, namely rotations in the {Pi,Px) plane. At 
this secondary fixed point, 6, the other "half" of the im- 
purity, acquires dynamics, rather than 7. The operators 
q, y, ij(3i^>l c ^ ec and ilPx^\ c ^ ec have scaling dimension 
3/2 and thus are irrelevant with RG eigenvalue -1/2. But 
of these, only the last two will give singular specific heat 
corrections (O(TlnT)). 

However, there is an important difference from the pri- 
mary 2CK fixed point discussed earlier. About this fixed 
point, there are two relevant operators, consistent with 
the symmetries of the model, namely i^Pi and ijPx (or, 
equivalently, i~/f3 c and ijP s ), both with dimension 1/2. 
These correspon d si mply to Ai and A_ before the change 
of variables (Eq(^)) leading to Eq(f4(i|); they can be gen- 
erated from nonzero q and y . The existence of two rel- 
evant operators is due to the fact that Pi and px trans- 
form the same way under the discrete symmetries; thus 
the artificial U(l) symmetry cannot be extended into an 
0(2) group (as was the case for the 2CK fixed point). 
As a result, this new fixed point is harder to find in the 
interchange symmetric case than the primary 2CK fixed 
point, as it requires the impurity-single-electron hopping 
term to be small and the two-electron- plus-impurity and 
the simple impurity hopping terms to be almost exactly 
equal at the Kondo scale. It should be noted, however, 
that the total number of relevant, dimension 1/2, oper- 
ators around this fixed point is again six, including the 
above mentioned ones, together with i^/d and the three 
spin SU{2) symmetry breaking operators, i^a s , ija es 
and ijPes, as discussed in Section IV. E and listed in Ta- 
ble |nj indicating that the fixed point symmetry is again 
that of the conventional 2CK model. This is supported 
by the fact that, just like in the case of the 2CK fixed 
point, there are four operators with scaling dimension 
3/2, albeit with completely different symmetries. 

Finally, some comments are needed on the the nature 
of this novel fixed point. We should first stress the ab- 
sence of Ai, the impurity-single-electron hopping term, 
which together with the Qo term (see Eq(||)) would form 
the conventional Kondo-like interaction term. Hence, 
the appearance of non-Fermi liquid behavior does not 
originate from the competition between the spin up and 
spin down electrons to form a channel-pseudo-spin singlet 
ground state with the impurity, but, rather, in the pres- 
ence of strong impurity-electron repulsion (Qo ~ 1/2), 
from the competition between bare impurity tunnelling 
and two-electron-plus-impurity tunnelling. Formally, 
there is an analogy with the conventional two channel 
flavor-anisotropic Kondo model,c3 with A± playing the 
role of the Kondo couplings of the two "flavor" channels, 
which can be seen in the left hand side of Eq(30) if $ c 
is substituted by <& s . When Ao ^ A2 one flavor channel 
couples more strongly to the impurity therefore screening 
it alone at low energy, which results in usual Fermi liquid 
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behavior. However, if Ao = A2 the flavor anisotropy dis- 
appears and the system flows to a non-Fermi liquid fixed 
point (provided Ai is zero). As a result, although this 
fixed point may be in the same universality class as the 
conventional 2CK model, the mechanism that brings it 
about is completely different physically. 



VI. ACCESSIBILITY OF THE TWO CHANNEL 
KONDO FIXED POINT AND CONCLUSIONS 



In the previous sections we have shown how the physi- 
cal operators in the tunnelling impurity problem behave 
near the two channel Kondo fixed point. In particular, 
we observed that a linear combination, A + = Ao + A2, 
of the bare impurity hopping, Ao, and the impurity-plus- 
two-electron hopping, A2, is relevant and drives the sys- 
tem away from this special critical point resulting in con- 
ventional Fermi liquid behavior at low temperatures, as 
for the usual one-channel Kondo system. If one could 
somehow tune A + , (or one other coupling such as the 
electronic hopping y) then one might be able to tune 
through the critical point and find the 2CK non-Fermi 
liquid behavior at low temperatures in the vicinity of a 
critical coupling. Unfortunately, such tuning over an ad- 
equate range is probably difficult to achieve. Thus one 
probably has to rely on the hope that a natural regime 
of couplings will lead to flow under renormalizatipn close 
to the 2CK fixed point. Vladar and ZawadowskH appear 
to suggest that this should be the case. Unfortunately, 
more complete analysis implies the converse, that only 
for fortuitous reasons would the impurity system — even 
without asymmetry between the sites — exhibit 2CK be- 
havior at low T. 

In this last section, we use the weak coupling analysis 
of Section III and Appendix A combined with the inter- 
mediate coupling analysis of Section IV and Appendix 
B, to find criteria for approaching close to the 2CK fixed 
point. 

Since we are interested in systems in which the Kondo 
temperature is much less than the bandwidth, the weak 
hopping behavior will control the relative strengths of 
couplings at the Kondo scale, at which the first of the 
hopping terms becomes of order the renormalized band- 
width. Operators which are irrelevant for weak hopping 
will flow away rapidly under renormalization, changing 
by finite amounts the remaining parameters, q, y and 
the {A^}, i = 0, 1, 2. For example, the complicated 
hopping-scattering term, w, which was discussed in Sec- 
tion IV and plays a role near the intermediate coupling 
fixed point, will be of order the hopping terms {A^} or 
smaller initially and flow away rapidly, modifying, among 
other terms, A + as in Eq(|4^), in the process. This will be 
the main role of such a term and we can incorporate its 
effects into a modified "bare" Ao and A2 . We thus start 
at an energy scale substantially below the bandwidth at 
which the important parameters for Q € [0, 1] are just 



q, y and the {A^} at this scale, the irrelevant operators 
having become small. The relevant eigenvalues for the 
hopping about the zero hopping fixed manifold will be 
universal. For small y they are given by Eq(|24|). 

If Q is initially small, i.e. q < 1/2, A+ will be sub- 
stantially larger than the other eigenvalues and thus a 
particular linear combination of the {A^} will grow fast. 
Unfortunately, as can be seen from Appendix A, this 
combination (Ao) is the wrong one to yield flow near 
the 2CK fixed point as it includes A+. Only if this com- 
bination, Ao, is initially very small relative to a power of 
another linear combination of the same {A^}, can behav- 
ior near the 2CK critical point be obtained; furthermore 
the criteria become more stringent the lower the Kondo 
temperature, as shown in Appendix A. 

Better prospects occur when Q as 1/2 (i.e. q small). 
Unfortunately, even if the one-electron plus impurity 
hopping term, Ai, were initially much bigger than Ao 
and A2, the purely electronic hopping term y — deter- 
mined basically by the spatial separation of the two im- 
purity sitesEJ — would combine with Ai to generate the 
wrong combination, A+, of Ao and A2, as in Eq(^3|). 
Thus, again, unless y is small the criteria from Appendix 
A are very strict, as could be anticipated from the y de- 
pendence of A± . 

The best prospects are thus for q and y both small so 
that the eigenvalues are all comparable. But this is just 
the condition for the analysis of Section IV and Appendix 
B via refermionization to be valid. To get near the 2CK 
fixed point, A + must rem ain small, thus we can study 
the RG equations Eq(B14) to leading order in A + and 
the linear combinations of y and q, y and q; 



dA 



K 



dl 
dq 
~dl 
d£ 
dl 



2 



1 - 



Air 



(44) 



A 2 K 

A 2 
8tt ; 



and, ignoring w from Eqd43) 



dA 



dl 2 A_I 



2yA K . 



(45) 



At the intermediate coupling 2CK fixed point, the con- 
vention we have chosen yields A^ = V4tt, so that y and 
q have the correct eigenvalues there as well as for weak 
coupling (Ak ~ 0). 

Integrating the above equations we find that the crite- 
rion to flow to the 2CK fixed point is, to leading order in 
A + and y 



A + -^ln 



1 



4-7T 



A 2 

K 



-in 



= 



(46) 



The parameters can be evaluated at any scale where the 
intermediate coupling RG equations are valid, i.e. for 
small A + , y and q. If these are small at the starting 
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scale, then their starting values can be used. Note that 

since Ak = \J A\ + (A — A 2 ) 2 , unless A and A 2 are 
almost exactly equal and opposite, Ak will be at least 
as big as A + initially so that |yln A^| needs to be small 
rather than just y, even in the best case of q = 1/2 — Q 
small (recall th at y is a linear combination of q and y 
given by Eq(B7)). 

We thus see that the flows near the 2CK fixed point, 
in particular the generation via Eq(ff3|) of the unique rel- 
evant operator, A + from other operators, yield stringent 
conditions for the accessibility of the non-Fermi liquid 
2CK behavior. If there is no symmetry between the two 
sites, then the presence of a second relevant operator (see 
Section IV) makes prospects even worse. 

This work strongly suggests that to observe two- 
channel-Kondo-like non-Fermi liquid behavior of an im- 
purity hopping between two sites in a metal one must 
either be able to tune some parameters over a substan- 



tial range, or be extremely hicky. This casts doubt on the 
interpretation of Ralph et an a of their narrow constric- 
tion tunnelling data. One possibility, although perhaps 
farfetched, is that these might be some kind of defects 
tunnelling in environments with higher symmetry, or at 
least approximate symmetry. In another paper, we will 
show how an impurity hopping among three sites with 
triangular symmetry can, without fine tuning, lead to a 
two channel Kondo behavior at low temperatures. 
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APPENDIX A: DETAILS OF WEAK-INTERMEDIATE COUPLING CROSSOVER 

In this Appendix we analyze the dependence of the crossover energy scale, Tk, and other quantities as a function 
of initial values of various coupling constants. For simplicity, we ignore the coupling to the even parts of the electron 
fields (and hence terms like u and w in Eq(p9|)). Since in the weak coupling RG flow equations (Eq(p3|)) of Ao, Ai and 
A2, q and y are essentially constant, these equations can be trivially integrated. The solutions for the two quantities 
of interest at the intermediate coupling region, i.e. A+ (I) and Ak (I) are 

A+ (0 = A (0 + A 2 (0 = A e x + l + A 2 e A - z (Al) 

and 



A K (I) = \J A\ (I) + Al (I) = J (A ie ^) +(A e^-A 2 e^) , (A2) 
with the weak coupling eigenvalues given by Eq(|24|) 

Ao = \ - 2q 2 

Note that the weak coupling eigenvalues for the relevant hopping terms should be universal.0 But since we have 
ignored higher order terms in y, (as well as renormalization of q from high energy scales which make it related to a 
phase shift rather than a coupling constant) the expressions for the eigenvalues in terms of the original y and q will 
not be exact. 

The couplings that appear in Eqs(Al) and (A2), Ao, Ai and A 2 , are linear combinations of the original Ao, Ai 
and A 2 : 

A = 2 ^-A_ ^ {(V^+^ + g) A + (v^T7-g) A 2 + yA 1 } (A3) 
Ai = = — = {y (A 2 - Ao) + gAi} 

A 2 = ^_{(v^T7-g) A + (>/?+7+g)A 9 -yA 1 }. 

The weak coupling flow equations in Eq(p3|) break down at the crossover scale Tk — r~ 1 e~ l , which is reached when 
the first of A + (1*) and Ak {I*) becomes of order unity. 
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At this scale, q and y, defined in Eq(B7) and appearing in Eq(flO|) are 



<1 



= vV + y 2 



(T K T C y° 

( A 



(A4) 



A, 



(T KTc y + (t kTc ) 



A_ 



In order to flow close to the 2CK fixed point, it is necessary that Ak (I*) = O (1) first, with A + (I*) <C 1. There 
are two regions of parameter space, in terms of Ao, Ai and A2, where this is realizable. The first is 





Ao 


Ao 







A+ 


Ai 






Ao 


A 2 








Ax 





< 1 



« 1, 



(A5) 



in which case the crossover scale is 



-K 



Ai 



(A6) 



The other region corresponds to the second square under the square-root sign in Eq([A2|) being of order unity. In 
this case the inequality of Eq(A5) is reversed, Aq and A2 have to be opposite in sign, (A0A2 < 0) and 







A+ 




A 2 












Ao 





- 1 



< 1. 



(A7) 



Then Tk is 



T, 



K 



Ao 



(A8) 



The opposite limit corresponds to A + (£*) becoming of order unity first. This will happen when Eqs(A7)-(A8) hold, 
but with Ao and A 2 having the same sign (A0A2 > 0). 

It can be seen from Eq(|A3j) that the above two conditions Eq(A5) and Eq(A7) are complicated functions of the 
initial parameters q, y and the A's. In the remainder of this Appendix we will analyze two limiting cases in terms q 
and y. 

The first limit is \q\ \y\. From Eq(A3) we see that for q > this cor resp onds to A, w A; for i — 0, 1 and 2, in 
which case the conditions in terms of Aj are read directly from Eqs([A^)-( A7). It should be noted that small initial 
Qo — i — <7o> i-e. weak coupling of the electrons to the impurity position, (Eq(||)), lies in this regime. 

In the second limiting case, whence |g| <C 



Ao 
A 2 
Aj 



2 

IA2-A0I 



(A + A 2 
(A + A 2 



sgn (y) Ai) 
sgn (y) Ai) 



(A9) 



In this region, EqflAq) becomes, for y > 0, approximately, 

|(A + A 2 ) + A 1 | < |A - A 2 
|(A + A 2 )-A 1 | < |A -A 2 



1 1+43/ 



(A10) 
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while Eq(|A7j), for y > 0, corresponds to 



|A 1 |>|A + A 2 | (All) 
|(A + A 2 ) - Ail 1 " 4 ^ « |(A + A 2 ) + Ail 1 ^ . 

In conclusion, for both limiting situations, as well as for the intermediate region, where \q\ w \y\, there is a finite, 
albeit small, region of initial parameters where the approach to the 2CK fixed point at low energies is indeed possible. 

Finally, we demonstrate that the weak coupling condition for flowing to the 2CK fixed point, A + (I*) small, evaluated 
at the Kondo temperature, Tk — r c ~ 1 e _! , matches with the one derived from the intermediate coupling side (Eq(|46[)) 
in the parame ter region where both are valid, i.e. small initial q, y, A + and Ak- To linear order in q, y and A + we 
have from Eq(|A2|), Ak(1) ~ e'^Aj? and hence 

T K nA 2 K T- x . (A12) 



The condition for criticality becomes from Eq(Al) 



—^4" + — = 0. (A13) 
{T K r c ) X+ (T K r c f- 

After expanding the factors T K X± to leading order in q and y and writing I* w — In (A^-) this yields 

( -^4r + — ^-jrl - 2^T^ ( -^4r - —^-A ln A k = o. (A14) 

\(T K r c ) Xo (T k t c ) Xo J \(T K r c ) Xo (T k t c ) " J 

Noting that the terms in the first parentheses is just A+ evaluated at Tk, to l owe st order in q and y, while the term 



multiplying the logarithm is just 2y at Tk, to lowest order in q and y (see Eq(A4)), we recover the condition Eq(|4f 
derived from the refermionization calculation, for the initial Ak small. 



APPENDIX B: DETAILED ANALYSIS OF INTERMEDIATE COUPLING POINT 



In this Appendix the details of the intermediate coupling point an aly sis are discussed. The starting point is the 
bosonization of the Hamiltonian in the Eqs(0-||), using Eq (|To|) and Eq([ll]), but inserting and keeping track of ordering 
factors which ensure proper anticommutation relations between different fields. The original fermions cif, c 2 |, c\y, 
c 2 j are bosonized in the following way (rather that just the simple form of Eq(|Io|)): 

(Bl) 



Cit = 


1 








C2T = 


1 








cil = 


1 






c 2i = 


1 


e 1 * 





where iVif, A^f etc. are the number operators for the associated fields: 

Nxi = I dxc{^ (x)cxi (x) = — I dec\ u cif e (B2) 



2tt . 

etc. with the total number of electron fields N tot = + A 2 | + N±i + N%i being always conserved. Note that since 
the number operators have integer eigenvalues the exponential factors in Eq(Bl) can only take the values plus and 
minus one. With these extra ordering factors the only change in the effective Hamiltonian of Eq(^) is the appearance 
of 

exp[i7r(7Vi T + A^ 2T )] (B3) 
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multiplying y and Ai. Since the combination + A^t is conserved in the absence of spin-flip processes, the 
exponential factor can be set to plus or minus unity and be disregarded. This can be seen more clearly if, using 
the number operators N c , N s , N ec and N es defined in an analogous way as the corresponding <I>-fields in Eq(|l9|) and 
Eq(20|) 1 we realize that +-/V21 = N ec + N es . Since N to t — 2N ec which is always conserved, the exponential factor of 



Eq(B3) will be conserved if N es is constant, i.e. if no spin-flip processes can appear. [Note that in the absence of spin 
flip processes, one can always choose a convention that different spin species commute rather than anti-commute.] 

After performing the unitary transformation using U = exp [— ^cr z ^c\ we get the transformed Hamiltonian Ti 
of Eq(p^). Next, following Emery and Kivelsonlia, we refermionize the exponentials of the Bose fields $ p with 



(j, = c,s, ec, es and also fermionizing er_ keeping track of ordering operators, as in Eq(Bl), to ensure proper anticom- 
mutation relations between different fermion fields {^^j and d: 

9 a = -l = e i *'e* KN ' (B4) 

V 2 7TT C 

^ _ 1 e i*c e i*(Nec+Ne.+Ns+N A ) 



d 



\J2tIT c 

jTr(N s +N es +N ec ) 



(T_e 



\J2ht c 
1 



V27TT C 

where is the number operator for the /i-fermions and = d^d. The Hamiltonian now becomes 

H = Ho + ^= { -A, (d + <fi) (*, - *t) + A _ (d + J) (* c - *t) } ( B5 ) 

+ ny {2d)d - 1) (tt c + *t) _ $t) 

- 2v q ^ c (2dU - 1) - (d - dt) (* c + * t) $t^ ec . 

Note that we have chosen the ordering operators in the refermionization, so that the ordering operators appearing in 
the bosonized Hamiltonian cancel. 

To simplify the Hamiltonian we introduce a set of eight Majorana fermions of Eq(|36|), so that the Hamiltonian 
takes the form of Eq(|3^). Their transformation properties under the various symmetries are listed in Table ||. Now, 
as it is clear from this Table, (3 S and (3 C have the same transformation properties under the discrete symmetries of 
the system. Therefore, we can make linear combinations of the two so that only one will couple to the 7 Majorana 
fermion. Defining 

A_/? c - Aifl. 

Pi = -r (B6) 



and 



Ai/3 C + A_/3 S 



y 



gA_ + yAi 



A 



K 



where Ak — J A\ + A 2 , as in Eq(|39|), we arrive to the Hamiltonian in the form of Eq(|40|) 

H = H + 2iua ec f3 ec + -J= {A K lPi + A+Sa c ) 

V27TT C 



+ AirqjSctcPx + 4iry~fSa c f3i + y/2irT c w5a c a ec (3 e 
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In the-RG the operator i"//3x is redundant; under renormalization it can be generated but can always be rotated 
away.a 

Before deriving the RG flow equations for the various operator amplitudes in Eq(f40|) close to the intermediate 
coupling fixed point, we should first locate the fixed point. As discussed in Section IV, formally it should lie at 
A a' = oo. This can be seen in the following way. At the fixed point, with A + = q = y = w = for convenience, the 
Ak term in Eq(|4C|) should have the same scaling dimension as Tio, i.e. unity. Defining the Fourier transform of the 
non-interacting (Ak — 0), imaginary time Green's functions for 7 and /3j as 

G 07 (iu) = ^- = -F.T. (T T [ 7 (r) 7 (0)]) (B8) 

iujg 1 

and 

G m (ia>) = -i S J^l = -F.t. {Tt [f3l ( T ) 0i (0 )]) o , (B9) 

respectively, where g 7 is related related to the normalization of 7, as discussed in Section IV, and F.T. stands for 
Fourier transform, we see that the corresponding interacting (Ak ^ 0) Green's functions are 

1 1 

Gy fa) = 7T2 = -7-2 7-T (BIO) 



G o 7 M " l^Gop, (%u) tg^u + V r ; 



and 



G Pl (tw) = ^ — = — 5^5— (Bll) 

This means that the naive scaling dimensions of 7 and /?/ (u —* 0) with A^- and <? 7 finite are 1/2 and 3/2, respectively, 
adding to 2. Thus some combination of Ak and g 7 must be infinite at the fixed point. But more useful than letting 
Ak — > 00 is to keep Ak finite and let <? 7 — > 0. A convenient form, which we will adopt, is to choose the normalization 
of 7 so that 

9-f + ^ = 1 (B12) 

This fixes the scaling dimension of 7 to be 1/2 at the fixed point, but also correctly gives the conventional scaling 
as Ak — > 0. This choice determines the part of RG flow equations for Ak and g 1 that does not involve any other 
parameters: 



dA 



K 



dl 

dfh. 
dl 



reseating 



reseating 



A 2 
-— ^0 



Thus, close to the weak coupling fixed point (Ak ~ 0, <? 7 =1), g y is marginal and A^- is relevant with dimension 1/2, 
while at the intermediate coupli ng fix ed point (A^- finite and g 7 ps 0), Ak is marginal and g-y irrelevant. In addition, 



the flow equation for in Eq(B13) determines the crossover scale, i.e. Ak (l*) — 0(1), where the renormalized 
values of the various amplitudes from the weak and intermediate coupling parts of the flows can be patched together. 

Next, we will derive the RG flow equations close to the 2CK fixed point. This involves four usual steps. We will 
use a simple frequency shell RG, taking the unperturbed Hamiltonian to be Hq with non-zero Ak and u; note that 
as usual, certain coefficients will depend on the RG scheme. First the modes with magnitude of their frequency in 
the infinitesimal range between t^ 1 (1 — dl) and r" 1 are integrated out. Then the frequency cutoff is rescaled to 
its initial value r" 1 and the fields are rescaled by their scaling dimension, determined by their two point correlation 
functions, from which the RG eigenvalues of the various operators can be determined. Third, perturbative corrections 
to the flow equations are included from integrating out the high energy modes in the four fermion terms, for small q, 
y and w. Thus, e.g. the contraction (or integration) of 7 and fti in the y-term contributes to the generation of A + . 
Finally, since at each step of this decimation process ij/3x will be generated, we have to rotate f3i and j3x so that it 
disappears. However, the corrections due to this rotation are negligible to the order at which we work. Following this 
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procedure we get the full set of RG equations describing the 2CK fixed point, and more generally for small A+, w, y 
and q, but with general A^, 



dA+ 

dl 
dA K 
dl 
dq 
~dl 

dy 
dl 
dw 

~dT 
du 
~dl 



= 1-A+ + 2yA K + 



wu 



2 

Ak 
2 



1 - 



K 



A 2 

A 2 
8tt ' 



A 

47T 

off 



%/2(l + u 2 ) 

+ 2yA+ + 0(y 2 A K ,A 

A+ _ 2 A+__A+ 
A^'A^A^ 



+ O (fA + ,w 2 A + ,wuA + ,A 3 + ) 



(B14) 



+ 



2 A+ _ 2 A + __A 
° {q ^ y ^ qy A K 



-- + 0(uA + ) 
wA + 

47T 



The constraint of Eq( |B12|) h as be en us ed; it has the effect of simplifying various messy terms that would appear in 
denominators from Eqs( |B10|) a nd ( Bll ) evaluated at |w| = 1/t c ££I The above flow equations are identical to the weak 
coupling flow equations Eq(|23D, to lowest order in y, q, Aj and qAi and yAi, with i = +, — , 1 or K . The higher order 
terms in the Aj's are non-universal and depend on the RG procedure. In the flow equations above, we have selected a 
RG scheme via Eq( B12| ), such that the scaling dimensions of the various operators are given correctly to leading order 
in the both weak and intermediate coupling limits. The u dependence of the w term in -^f- can, however, depend on 
the RG scheme, or, equivalently, the definition of u. 

If w, A + and y are all small at some scale, then the wA + and yA + terms can be ignored. The RG equations can 
then be integrated to yield the condition for lying on the critical manifold: 



A i 



K 



V2(l + u 2 ) 



2yA 

A 2 

4-7T 



In 



A 2 



4tt 



= 



(B15) 



to linear order in w, y and A + . In general, the parameters above should be evaluated in the regi me of the crossover to 
intermediate coupling. [In the approximation of ignoring the neglected terms, the condition Eq(B15) can be used at 
any scale.] The role of w in the above equation is a typical example of the way an irrelevant operator can renormalize 
the value of a relevant one before decaying away. In this context, one should remember that the values of both Ak 
and A + should include corrections from the more strongly irrelevant operators that exist at higher energy scales but 
become negligibly small as we approach the fixed point, or as the flow moves away from weak coupling. 

The effects of various symmetry breaking terms are discussed in Section IV. D. The symmetries of the important 
operators are summarized in Table O. If spin-flipping is allowed, it was shown in Section IV. D that the relevant 
processes involve joint electron-impurity hops which carry electronic spin. Their representation in terms of the bare 
fields involves ordering operators: 



-A' es i{a + e^-a^e-^ 



A s 



[u+e — cr_e 



— i$ ea \ 



(B16) 



)} 



but after refermionizing, these disappear with the choices made in Eq( |B4|) , yielding, in the Majorana representation, 

Hsf = 2V2~¥T c i (A es Sa es + A' es 6(3 es + A s Sa s ) . (B17) 



APPENDIX C: LEADING IRRELEVANT OPERATORS - COMPARISON TO CONFORMAL FIELD 

THEORY RESULTS 

In this-gaper we have implemented an RG approach based on a bosonization method first developed by Emery and 
Kivelsor£3 which is inherently anisotropic in the channel pseudo-spin. We thus found a 2CK-like fixed point with an 
0(2) symmetry, higher than the Z2 interchange symmetry of the original impurity problem. Despite the fact that this 
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fixed point had lower symmetry than the SU{2) symmetric one of the 2CK model solved by conformal field theory 
methods, we will show in this Appendix that all relevant, marginal and leading irrelevant operators expected from 
conformal field theory can be recovered from this model.113 

To be specific, in Section IV. D the six dimension 1/2 relevant operators of the 2CK fixed point were analyzed (see 
Table p|) . Also, the only marginal, dimension 1, operator was found to be the conventional potential scattering term, 
u^lc^ec, yielding ia ec f3 ec . 

In addition, there are four leading irrelevant, dimension 3/2, operators consistent with the symmetries of the 
problem we considered. The first, ^ySctcPx, corresponding to q, does not break the 0(2) symmetry of the fixed point. 
This is essentially the, leading irrelevant operator discussed in reference [h]. In fact, there is no other operator with 
the same dimensionJHj invariant under 0(2) and the other discrete symmetries of the system, that does not involve 
even operators (^ec, & e s)i which naively decouple from the impurity. This can be checked directly from Table |, by 
taking all possible combinations of Majorana fermions in the channel pseudo-spin sector. Since the 2CK fixed point 
can be correctly described in terms of the free Majorana fermions we conclude that there is no other operator missing 
from our approach. This fact suggests that the operators that would distinguish the 0(2) from the SU(2) symmetric 
fixad points have dimension higher than 3/2 and are thus more strongly irrelevant, in contrast with what Affleck et 
a&A have claimed that a term breaking the SU(2) down to an 0(2) symmetry has dimension 3/2.EZI 

The second leading irrelevant operator, b~a s a es (3 es , can be seen from Table § to be the only other operator with 
dimension 3/2 which respects the 0(2) as well as the discrete, symmetries of the system. Although it explicitly 
involves the Majorana fermions of the spin sector, it can be seen that this operator also satisfies the full spin SU (2) 
symmetry by the fact that it is generated from a strictly spin SU(2) invariant term, Ai. Starting from Eq(||), one 
can first bosonize the fermion operators c\ a directly into exponentials of combinations of $ e c, &es, 3>c and $ s . Then 
the Ai operator becomes 



Ai 

AlTT a 



e -£(*«e+*«-*o-*,) e 4(*ec+*e,+*o+*,) (£1) 



(multiplied by expi7r {N es + iV ec ), the appropriate ordering operator discussed in Appendix B). Obviously, the leading 
operator of this term is the one appearing in Eq(|l3[). However there are subleading corrections, coming from the 
operator product expansion of the various exponentials. One such term is 

i (cr+e**" - (T-e-**-) ^sin$„ (C2) 

with a cutoff dependent coefficient proportional to Ai. This is allowed by the full symmetries of the 2CK problem 
and will therefore be generated in the way outlined above or, equivalently, from the high energy scattering processes. 
At weak coupling, this term, with RG eigenvalue —1/2 — 2q 2 , is always irrelevant. After performing the unitary 
transformation of Eq(|28|) and refermionizing, it becomes 5a s a es {3 es , with scaling dimension 3/2 close to the interme- 
diate coupling fixed point. This operator contributes to the singular part of the impurity specific heat but not to the 
susceptibility of the conventional two channel, spin 1/2, Kondo model, thus leading to a non-universal Wilson ratio 
(The existence of two irrelevant operators at the 2CK fixed-point was discussed by Affleck and Ludwig in Appendix 
D of reference 12). Other approaches to the 2CK problemEj'liSEa assumed implicitly some higher symmetryc£l at the 



fixed point and therefore do not contain this second leading irrelevant operator (see reference [Tl] for more details). 

The remaining two dimension 3/2 operators, j5a c /3i and 5a c a ec /3 ec , do not respect the 0(2) symmetry, but exist 
on the critical manifold with the Z2 interchange sym metry , i.e. having set the relevant operator, corresponding to a 



linear combination of A + , w and y, to zero — see Eq(B15). They were discussed in Section IV. C and correspond to 



couplings y and w. Interestingly, the first does not enter in the leading impurity specific heat correction of 0(T In T), 
since it couples to a dimension 2 irrelevant operator, g 7 , leading only to a T 3 lnT correction to the impurity specific 
heat. 

More generally, in the case where Z2 interchange symmetry is also broken, there are three more dimension 3/2 
irrelevant operators, namely r ySa ec /3 ec , iSQgjr and iS^-. However, in this case two relevant operators must be tuned 
to be on the critical manifold. 

Finally, for completeness, we summarize the leading irrelevant operators for the additional non-Fermi liquid fixed 
point discussed in Section V. Apart from q and y, the operators 7/3/a e c/3 ec and ~/Px a ecPec also have scaling dimension 
3/2. Again, only the last two contribute to 0(T In T) to the impurity specific heat. As in the above case, in the absence 
of Z2 symmetry, three additional dimension 3/2 operators will be present, namely *7^S jd&ecPec and ja c f3 c f3 s . 
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Fermion 
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Pes 


+ 
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TABLE I. Transformation properties of Majorana fermions with respect to the discrete symmetries of the system. Electron 
number conservation implies that the Bose field can be shifted by an arbitrary amount, which means that only the 
combination a ec /3 ec can appear. In the presence of conservation of the ^-component of the spin, "I> es can also be shifted by an 
arbitrary constant. If z-spin were only conserved mod 2, then just a discrete symmetry would remain; this is listed in the last 
column of the table as a simple way to distinguish the different symmetries of the e operators. 



Relevant Operators 


Site Interchange (Z2) 


Time Reversal 


Spin Reversal 


Spin Conservation mod 2 


iSa c 
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iSPx 
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i5a a 
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iSa e3 




+ 






iSPes 






+ 





TABLE II. The six independent relevant operators at 2CK fixed point and their transformation properties (for last column, 
see note in Table I). Note that iSPs and iSp c are not listed here but rather the combination iSPx which has the same symmetries 
as the other two. (See discussion in Section IV. D). 
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